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PREI ACE 


Date, selection and reduction procedures are described by 
which scalar and vector magnetic anomaly maps are constructed. 

The scalar and vertical magnetic anomalies are believed to be 
generated mainly in the earth's crust. The horizontal anomalies 
are believed to be mainly due to short-period field-aligned cur- 
rents. The correlation of scalar magnetic anomalies with known 
oceanic stricture is remarkable — magnetic highs are associated 
with oceanic ridges and magnetic lows with abyssal plains. The 
correlation teeween anomalies and continental geology is not as 
clear. In East Antarctica, magnetic lows associated with the 
Ross Embayment and the Amery Ice Shelf are consistent with the 
hypothesis (Hayes and Davey, 1975) that these regions are failed 
continental rifts. The magnetic low over the Gamburtsev Moun- 
tains, as distinct froir the highs over Wilkes Land and Enderby 
Land, may imply a tectonic history of the Gamburtsev Mountains 
that is different than that of the latter two surrounding regions 
(Drewry, 1975). The Wilkes Land magnetic high is consistent with 
the hypothesis (Veevers, 1982) that this region is the site of 
cool convergence in the mantle. In West Antarctica, Dalziel and 
Elliot (1982) have postulated the existence of five microplates. 


Separate anomaly features are associated with each of their micro- 
plates, supporting their hypothesis. Finally, the apparent 
mirroring of anomaly features on the other Gondwanaland conti- 
nents implies that the sources of the anomalies have been stable 
since before the rifting occurred. 
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I. INTRODUCTION 

It is surprising and fortunate that magnetic fields generated 
by crustal sources can be isolated in satellite magnetic field 
data from those generated in the core and those external to the 
earth. It is surprising because crustal fields (0-25 nT) are so 
much smaller than the core field (30,000-60,000 nT) and external 
fields (0-2000 nT) , and fortunate because geologic and tectonic 
features in the deep crust can result in long-wavelength magnetic 
anomalies (Pakiser and Zietz, 1965; Zietz et al, 1966: Hall. 

1974; Krut ikhovskaya and Pashkevich, 1977). Long -wave length mag- 
netic anomalies (> 200 km) were first mapped globally from data 
acquired by the POGO (Polar Orbiting Geophysical Observatories) 
satellites (Regan et al, 1975) and have compared favorably with 
upward-continued regional aeromagnetic surveys (Langel et al, 
1980). The chief advantage of the MAGSAT data over the POGO data 
lies in the ability to study vector fields with MAGSAT and in the 
greater resolution MAGSAT' s generally lower orbit affords 
(perigee 352 km, apogee 561 km). A disadvantage for Antarctic 
studies, outweighed by these advantages, is that MAGSAT' s orbit 
is inclined at 7° off the earth's axis of rotation resulting in a 
data gap around the Pole in MAGSAT derived maps. 

Attempts have been made to isolate anomalies both in the 
scalar field (Langel et al, 1982a; Coles et al, 1982; Ritzwoller 
and Bentley, 1982 (Appendix I)) and in the vector field (Langel 



et al, 1982b; Coles et al, 1982). Sailor et al (1982) confirm 
the general reliability of such attempts and conclude that resolu- 
tion is possible down to a 250 km spatial-wavelength in mid- 
latitudes; this being larger in higher latitude regions due to a 
smaller signal-to-noise ratio. The reduction in this ratio is 
due to the amplification of external magnetic fields in auroral 
regions. (Appendix IV describes the structure of current systems 
external to the earth and their implications to MAGSAT studies.) 
Magnetosphcric ring-currents in low latitudes generate a temporal- 
ly rather stable long-spatial wavelength signature (Langel and 
Sweeney, 1971) that is easily filtered. However, currents follow- 
ing magnetic field lines in auroral regions have temporal and 
spatial spectra covering quite a broad band of frequencies, thus 
making standard frequency-based filtering difficult. Therefore, 
only passes occurring during low field-aligned current periods 
can be used as data in high latitudes seriously reducing the size 
of the data set over Antarctica and in the Arctic (Coles et al, 
1982). The map of Ritzwoller and Bentley (1982) over Antarctica 
demonstrated the signal detection problem by displaying serious 
non-physical radial striping. 

Progress in modeling and interpreting satellite magnetic 
anomalies is short and generally qualitative. On a global scale, 
anomalies appear to be associated with such large structures as 
continental shields and platforms, subduction zones (positive), 



oceanic ridges (positive), basins (negative), and abyssal plains 
(negative) (signs are for reduced-to-pole maps) and appear to be 
bounded by such "linear" features as sutures, rifts, folded moun- 
tains, and age province boundaries (Frey, 1982a). More detailed 
preliminary regional studies have been performed by Frey (1982b) 
for Asia, by Hastings (1982) for Africa, by Hinze et al (1982) 
for South America, and by Ritzwoller and Bentley (1982) for Ant- 
arctica. 

A geologic interpretation of the maps must only be conducted 
in light of the probable mineralogy of the lower crust and upper 
mantle. Wasilewski et al (1979) argue that the mantle is proba- 
bly nonmagnetic so that if the Curie isotherm is below the crust 
the lower magnetic boundary is the Moho. Moreover, Wasilewski 
and Mayhew (1982) conclude that for some tectonic settings, at 
least, the lower crust may be the most magnetic crustal layer and 
that magnetization values for lower crustal xenoliths (specifical- 
ly metabasic rocks of the granulite facies) have values consist- 
ent with those inferred from models of long-wavelength anomalies. 
Wasilewski and Fountain (1992) corroborate these findings with a 
study of the Ivrea Zone in Northern Italy, wh mafic granulite 
facies rocks are the only magnetic lithology present, are thick 
and laterally continuous, and thus provide a good candidate for a 
deep-crustal source of long-wavelength magnetic anomalies. 


Along with regional mineralization variations, it is suspect- 
ed that variations in the depth to the Curie isotherm (if above 
the Moho) will contribute to long-wavelength anomalies (Mayhew, 
1982). Thus, the main sources cf long-wavelength magnetic anoma- 
lies in continental regions ai expected to reside above the Moho 
and tne Curie isotherm but principally in the lower crustal layer. 
Consequently, since the causes of variations in heat flow and 
regional mineralization are many, ambiguity is injected into the 
interpretation of satellite magnetic anomaly maps. However, 
everything else being equal, heat flow and regional magnetic 
anomalies are inversely related, wherecs the thichness o. the 
magnetized crust and the regional magnetic anomalies are d. . -y 
related. Therefore, continen:al highs indicate either a th. ~ 
crust with a deep Curie isotherm (low heat flow), or exceptional- 
ly high magnetic susceptibility in the lower crust, or both. On 
the other haM, continental lows imply some combination of a thin 
crust, a shallow Curie isotherm (high heat flow), and low suscep- 
tibilities in the lower crust. An additional degree of ambiguity 
in the determination of the natuie of the sourc; of an anomaly is 
also added by the large areal extent of the potential source re- 
gion. It remains uncertain to what degree continental remanent 
magnetization may affect long-vave length magnetic anomalies. 
Wasilewski and Padovani (1980) argue that the effect should be 
minimal, but Galliher and Mayhew (1982) were not able to demon- 
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strate that fact from an equivalent-dipole analysis of data over 
the United States. 

In the oceanic crust, the highest susceptibility region is 
in Layers 2 and 3, much nearer to the surface than in the conti- 
nental crust, and generally well above the Curie isotherm, except 
right at spreading centers. Thus, oceanic magnetic anomalies 
should reflect regional crustal thicknesses and susceptibility 
differences, but generally will not reflect heat flow except in 
extreme circumstances. Therefore, oceanic basins with a thin 
magnetized crust are expected to display a negative anomaly rela- 
tive to continental regions where the crust is generally thicker; 
but oceanic ridges can yield magnetic highs due to remanent mag- 
ne. zation even at satellite elevations, as has been shown by 
LaBreque (personal communication, 1983). Topographic variations 
near ridges, assuming induced magnetization alone, cannot account 
for the satellite magnetic highs associated with them. 
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II. DATA REDUCTION 

Ritzwoller and Bentley (1982) described a method by which a 
preliminary map of the «.rustal scalar magnetic field was produced. 
The method was designed to filter non-crustal magnetic fields 
using scalar magnetic field data alone. However, since a great 
deal of information concerning the strength and direction of 
field-aligned currents is available in the vector magnetic-field 
data, the optimal data-reduction procedure must consider the 
vector data. With this in mind we have followed the data selec- 
tion and reduction procedure described below in constructing the 
magnetic anomaly maps found m Figures 5-8. 

After the NASA Inves t igator-B data tape had been translated 
into Harris-computer compatible code by the program IBMTR (all 
computer programs referred to in the text can be found in 
Appendix V) and reformatted by REFMT, the data were in a shape to 
be reduced and selected. If a pass over Antarctica took place 
during a relatively magnetical ly-quiet time period (planetary 
magnetic activity index, Kp, less than or equal to 1” for 6 
hours), the degree and order 13 spherical harmonic core-field 
model MGST(4/8l) created by Langel et al (1980) and briefly des- 
cribed in Appendix III was subtracted from the scalar magnetic 
data calculated from the observed vector magnetic data. (The 
scalar magnetometer onboard the satellite malfunctioned, making 
it necessary to calculate the scalar field from the vector data.) 
Of the approximately 2400 passes over Antarctica between Novem- 
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ber 1, 1979 and April 1, 1980, 212 met this selection criterion; 
the remainder of the passes were discarded. To these passes a 
quadratic polynomial was least-squares fitted to the resulting 
spatial-series in an attempt to filter the effects of magneto- 
spheric ring-currents, core field model bias, and other errors in 
measurement. (For more details see Langel et al, 1982a). 

The vector and scalar magnetic anomalies were, thus, cal- 
culated in the following way. Let X£, Y£, and Z£ be spatial 
sequences of observed vector magnetic anomaly data in the right- 
handed rectilinear coordinate system (x,y,z) in which the posi- 
tive x-direction is North, y is East, and the z-direction is 
positive down along a radial vector connecting the satellite with 
the center of the earth. Then the vector anomalies in each direc- 
tion are the sequences AX£, AY£, and AZ£ calculated as follows: 



AX£ 

- X£ 

- XMOD£ 

- fi 


(la) 


AYi 

= Yi 

- YM0D£ 

- fi' 


(lb) 


AZ£ 

3 Zi 

- ZMOD£ 

- f£" 


(lc) 

The 

sequences 

XMOD£ , 

YMOD£, and 

ZMOD£ represent 

the vector core- 

field model values in 

the x,y , 

and z directions. 

respect ive ly , 

and 

the 

sequences f £ , 

f £ ' . and 

f £ 11 are quadratic 

polynomial 


sequences each least-squares fitted (Levington, 1969) to the 
first two terms of the right-hand-si le of (la), (lb), and (lc), 
respectively. The quadratic polynomials were fit to the se- 
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quences in the hope of removing the long-wavelength magnet ospher- 
ic ring-current effect described by Langel and Sweeney (1971). 
Magnetoapheric fields tend to impart a random "level" offset and 
bias to each profile that renders the data set inconsistent. It 
has been shown (Regan et al, 1975) that subtracting a best fit 
linear, quadratic, or first zonal harmonic fit over a limited 
latitude range makes the data set internally much more consistent. 

The scalar magnetic anomaly sequence, AB{, was calculated as 
follows : 

ABf */x£ 2 + Y£ 2 + Z£ 2 -V / XMOD i 2 + YMOD£ 2 + ZMODj 2 - g£ (2) 

where gj is the quadratic fit to the first two terms on the 
right-hand-side of (2). Figures la-d are graphs of the anomalies 
created by this process. This procedure was carried out by the 
program ANTAP. The data were then gridded by the program MRGRD. 

At this stage in the processing the majority of the passes 
still showed field-aligned current effects in the vector data. 
Figure 1 demonstrates what these effects "look" like in the 
vector data. Since field-aligned currents, the assumed cause of 
these high amplitude disturbances, are nearly normal (within 20°) 
to the satellite's path within the region of high field-aligned 
current effects, we see from Maxwell's first law: 
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Figure 1: Example of data from a pass (Pass 51) believed to be 

relatively uncffected by field-aligned currents. Graphs from top 
to bottom: scalar anomalies, AB; vertical anomalies, AZ; radial 

anomalies, AX; and tangential anomalies, AY. The AB graph explic- 
itly shows the polynomial fit to the data sequence corrrecting 
the "level offset". For simplicity the vector graphs leave this 
out . 


I 
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Vx AB' « J 


(3) 


where J represents the time-varying field-aligned current and AB' 
the magnitude of the magnetic field induced by this current, 
that : 


Vx VB' * J z z (4) 

where Jz is the component of J in the z-direction and z is the 
unit vector in the z-direction, so that if AB' * (AB'x, AB'y, 
AB'z) then: 

- 7 - AB'y - — - AB'x * yj (5) 

2 x 2 y z 

That is, the majority of the magnetic field arising from field- 
aligned currents is to be found in the horizontal magnetic field. 
This is observed in Figure 1. 

A couple of characteristics of the vector anomalies are 
worth noting. First, the location of the peaks in AX, AY, and 
AZ (here and afterwards subscripts will be suppressed) are highly 
correlated; notably, high-amplitude disturbances in AZ (> 25 nT) 
are always at points at which AX and AY are also disturbed but 
with the latter fields having higher amplitudes (50-300 nT). 

With respect to the horizontal field, a peak in AX often corres- 
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ponds to a peak and a trough pair in AY. This effect can be 
understood if we view the source of the disturbance as a line 
current nearly normal to the flight path. Then AB' , resulting 
from Jz, will reside in the horizontal plane as shown in Figure 2 
in which the satellite is at point A. As the satellite passes 
through quadrant III of the current's coordinate system, AX and 
Ay will show effects due to Ab'. Specifically, AX and AY will 
both show negative anomalies due to AB 1 . However, in quadrant 
IV, AX will show a negative anomaly but AY will show a positive 
anomaly due to AB'. Table I summarizes the behavior of AX and AY 
resulting from crossing a given quadrant of the current's coordi- 
nate system. Thus, the disturbances in AX and AY can be under- 
stood as arising from multiple line-current sources, consistent 
with the belief that field-aligned currents are generating the 
horizontal disturbances. 

Second, it is worth noting the relationship between the 
scalar anomalies and the vector anomalies from which they are 
calculated. Since positive AZ anomalies are directed toward the 
earth's surface, the expected positive correlation between the 
scalar anomaly and the vertical component of the vector anomaly 
would be expressed in Figures 1, 3, and 4 as negative correla- 
tions between AB and AZ. (In Figure 8, signs have been reversed 
to bring out the relationship with AB.) Thus, due to the nature 
of the coordinate system (x,y,z) a highly magnetized region of 
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Figure 2: Schematic representation of the magnetic field, AB 1 , 

generated by a field-aligned current, Jz. The satellite, A, 
follows the x-axis and crosses quadrants III and IV in the 
current-centered coordinate system. The effect that A B' has on 
AX and AY is summarized in Table I. 
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TABLE 1 


Effect of a Field-Aligned Current on the Horizontal 
Anomalies in the Current's Four Quadrants 


Quadrant 

I 

II 

III 

IV 

AX 

Positive 

Positive 

Negative 

Negative 

AY 

Positive 

Negative 

Negative 

Positive 
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Figure 3: Example of data from a pass (Pass Zl 7 ) shoving signs 

of a long-wavelength field-aligned current perturbation apparent 
on the scalar as well as the vector d»ca. Short -wavelength varia- 
tions apparent in the vector data do not appear on AB. This pass 
was discarded from the data set. 


* 
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Figure 4: Example of data from a pass (Pass 1887) shoving high 

amplitude short-period field-aligned current perturbations in the 
vector anomalies which do not appear on AB. However, long- 
wavelength perturbations due to field-aligned currents are 
believed to exist in AB since this pass is poorly correlated 
with nearby passes. This pass was also discarded from the data 
set and illustrates the necessity of the use of vector anomaly 
data to identify the effects of auroral currents. 
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the crust would tend to increase AB and decrease AZ. Peculiarly, 
though in Figure 1 AZ is roughly anti-correlated with AB, it also 
shows larger high-frequency components and is generally higher in 
amplitude. Furthermore, the scalar anomalies in every pass do 
not show the high frequency variations which are observed in the 
vector anomalies. Data processing has been carefully checked and 
no unwanted averaging or low pass filtering has been performed. 
Furthermore, the effect is not observed at lower latitudes so 
that magnetometer error is unlikely. Thus the effect appears to 
have a physical source, presumably the field-aligned sheets (see 
Appendix IV). To a first approximation these may be viewed as 
carrying the same current d nsity, and as being of large 
horiicn < l extent. Since the current is upward in one sheet and 
downward in the other, and since the magnetic field caused by a 
current sheet at a distance small with respect to the horizontal 
extent of the sheet depends only upon the current density and not 
upon that distance, the net effect observed in passing through 
the paired sheets would be a rotation of the ambient magnetic 
field without a change in its magnitude. 

Interestingly, longer-period variations in AZ are observed 
in AB but the short-period variations superimposed on these are 
not observed (Figures 3a-d). Apparently either the fields 


produced by the two current sheets are not the sa. 11 * at longer 
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periods, or the longer-period effects result from currents 
elsewhere in the ionosphere or magnetosphere. 

These examples show that studies performed with a scalar 
magnetometer alone will not be able to determine if moderate 
field-aligned currents are affecting a pass. Thus, the scalar 
anomaly map of Ritzwoller and Bentley (1982) is probably affected 
by field-aligned currents whose signatures were indistinguishable 
from crustal signatures in this scalar study. However, if by use 
of a vector magnetometer one can determine that high -amplitude 
rapid variations exist in AZ, then field-aligned currents are 
affecting the data and the pass can be deleted even though the 
scalar profile looks quite calm (Figures 4a-d). The pass repre- 
sented in Figures 4a-d would provide undesirable scalar data. If 
the anti-correlation between AB and long-period variations in AZ 
is quite high, then the effect of field-aligned-currents can be 


taken as minimal. 
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With this in mind we have further required each pass to 
satisfy the following selection criteria to be accepted as an 
"unaffected" pass: 

a. cne maximum amplitude of the vertical vector 
anomaly, AZ, is less than 25 nT 

b. the maximum amplicude of the scalar anomaly, 

AB, is less than 20 nT 

c. AB and -AZ are highly correlated. 

The data selected in this manner appear to strike the bal- 
ance between the procession of desirable data density (approxi- 
mately one pass per 2° on the average) and the minimization of 
field-aligned current effects. The scalar and vector magnetic 
anomaly data from the 88 passes satisfying these criteria were 
separately a/eraged in square bins measuring 330 km on a side and 
standard deviations were calculated. This was performed by the 
programs BIN and RDBIN. The data density for the 24 x 24 array 
centered on the geographical South Pole is found in Table I, 
Appendix VI. Tracks cf these 88 passes can be seen in Figure 10 
and pass numbers, dates, times, magnetic activity indices, and 
average elevations for each pass are listed in Table II, Appen- 
dix VI. Since NASA's pass numbering system begins each pass in 
the middle of Antarctica, each flight over Antarctica contains 
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Figure 5: MAGSAT total-field magnetic anomaly map over Antarctica. 

Units in nT. Average elevation 470 km. Capital letters indicate the 
approximate location of: A - Bellingshausen Abyssal Plain; B - Weddell 
Abyssal Plain; C - Enderby Abyssal Plain; D - Wilkes Abyssal Plain; 

E - Maud Rise; F - Kerguelen Plateau; G - juncture of the Mid-Indian 
Ocean Ridge and the East Pacific Ridge; H - South Sandwich Islands; 

I - Antarctic Peninsula; J - Ellsworth Mountains; K - Queen Maud Land; 

L - Enderby Land; M - Amery Ice Shelf/Lambert Glacier; N - Wilkes Land; 
0 - Transantarct ic Mountains; P - Ross Sea embayment ; Q - Marie Byrd 
Land; R - Thurston Island; S - Gamburtsev Mountains; T - Weddell Sea 
embayment . 


<♦ 
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Figure 6: Radial, AX, anomaly map. Positive field directed 

radially away from the South Pole at each point. Average 
elevation. 470 km. Average standard deviation: 4.3 nT. 

Contour interval: 10 nT. 


1 
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Figure 8: Vertical, AZ, anomaly map. Signs of anomalies 

r°*,ersed from the profile data so that positive is directed away 
from the earth's surface. Average elevation: 470 km. Average 

standard deviation: 2.0 nT. Contour interval: 2 nT. 
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data from two of NASA's passes. We have adopted the following 
convention in referring to pass numbers: our pass n contains data 
from NASA's south-going pass n and north-going pass n + 1. How- 
ever, when referring to other characteristics of our pass n over 
Antarctica, such as dates, times, etr., we have chosen these 
values from NASA's pass n + 1 si:»ce these values were determined 
when the satellite was at its southern node. 

In constructing Figures 5-8, ■> f any AB, AX, AY, and AZ values 
within a bin departed, respectively, by more than 4, 10, 10, and 
5 nT from the average for that bin, these data points were reject- 
ed and the averages and standard deviations recomputed. Between 
10 and 152 of the data were rejected on this basis. The average 
AB, AX, AY, and AZ anomaly values for each bin in the 24 x 24 
array can be seen in Tables III-VI, Appendix VI. 

If data are nearly normally distributed within each bin, 
magnetic values will cluster around the bin-average which is 
hoped to be close to the magnitude of the crustal field. Field 
aligned currents will create data values well away from the bin- 
average, though to the degree these effects are random around the 
bin-mean, averaging will diminish their effect by cancellation. 

A map of the scalar-map bin-standard-deviations, calculated be- 
fore discarding data based on the 4 nT discards 1-criterion , shows 
that in 3 of the 4 quadrants the anomalies of Figure 5 are based 
on data with standard deviations less than 3 nT. However, in the 
grid-southeast quadrant, the quadrant nearest to the geomagnetic 
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South Pole, standard deviations are significantly higher, ranging 
up to 6 nT. Thus, the features in the grid-southeast quadrant of 
Figure 5 have to have larger "mental" error-bars assigned to them 
than the features in the remaining quadrants. As will be discuss- 
ed in section III, possible perturbations due to the horizontal 
anomaly can be seen in Figure 13. The continental portion of the 
grid-southeast quadrant is relatively unaffected by the horizon- 
tal component of the magnetic anomalies. Thus, the large vari- 
ances in this quadrant must be due to variations in AZ. However, 
discarding data more than 4 nT from the bin-average acts to de- 
lete "affected" data points so that the recalculated bin-average 
should be closer to the value of the crustal fieid. This reduces 
the size of errors in all quadrants appreciably. Non-normal dis- 
tributions of data within a bin (e.g. bimodal distributions) will 
make this technique sub-optimal, however visual scrutiny of the 
data indicates that most values do lie around the mean with a 
smaller percentage of outliers on either side. 

Finally, the bin-values were hand-contoured yielding the 
maps in Figures 5-8, in which the average bin elevations for each 
map is approximately 470 km and the average bin standard varia- 
tions are 1.5, 4.1, 4.3, and 2.0 nT, respectively. The AZ map in 
Figure 8 contains a sign reversal so that the correlation with 
the scalar anomaly map will be clearer. The subtraction of the 
polynomial makes the map-average the zero-level on each map. 
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Thi- zero-level is arbitrary but reasonable. Since the map is 
comprised of approximately half continental and half oceanic 
regions, the former on the average higher in magnetic amplitude 
than the latter, this choice of zero-level approximately bisects 
the interval between oceanic basin lows and thick continental 
highs. This is a desirable location for the zero-level. Notice 
that in Figure 5 the zero-contour approximately corresponds to 
the edge of the continental shelf between the Wilkes Land magne- 
tic high and the Wilkes Abyssal Plain low, as we would desire. 

The anomalies on these maps have not been reduced-to-the- 
pole nor have data been continued to a single elevation before 
processing. However, since all of Antarctica is above 60°S geo- 
magnetic latitude, reduct ion-to-the-pole will have minimal effect. 
Furthermore, studies performed elsewhere (R. Sailor, personal 
communication, 1982) indicate that maps created in the w^y de- 
scribed above are "qualitatively and quantitatively similar" to 
those continued to the average-bin elevation from the same data 
set, using the equivalent source technique of rtayhew (1979). 

That is, effects due to data points that are higher than average 
approximately balance the effects due to points lower than aver- 
age so that their net effect on the average is small. Since 
fewer passes meet the selection criteria at high latitudes the 
data set is less dense so that cancellation will be less perfect 
than in lower latitudes. A map consisting of average-bin eleva- 
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tions shows considerable variability, but the elevation varia- 
tions appear to be unrelated to variations appearing on the 
scalar anomaly map in Figure 5. Furthermore, the data set has 
been partitioned into two subsets, -,»ie containing data from 
passes with an average elevation less than 475 km and the other 
containing data with average elevations greater than 475 km. 

Each subset contains 44 passes. The maps are highly correlated 
with each other and with Figure 5, though, as expected, the lower 
elevation map shows anomalies with somewhat larger amplitudes 
than the anomalies in Figure 5. These tests appear to indicate 
that elevation variations in the data have not generated an 
appreciable error in Figure 5. 
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III. ORIGIN OF THE ANOMALIES 
A. Scalar Anomalies 

The reduction of the data consisted of modeling magnetic 
fields generated in the core and magnetosphere and discarding 
passes which showed signs of severe field-aligned current effects 
registered in the vector anomalies. Any remaining field-aligned 
current effects random within the five month data window are 
hoped to be stacked out in creating the 3° square bins. How 
successful have we been at eliminating fields from each of these 
sources in the scalar anomaly map? 

First, considering the accuracy of the core-field model, 
Carle and Harrison (1982) have argued that the field remaining 
after removal of the spherical harmonic core-model will have a 
significant long-wavelength component. The majority of the core 
field rtmaining after removal of the model field is of very long 
spatial-wavelengths and is removed by the polynomial fit to the 
data. It is unclear what percentage of the remaining anomalies 
is due to the core, though it is reasonable to believe that a 
fraction of very long-wavelength anomalies (> 3000 km) has its 
origin in the core. We believe that this fraction is insignifi- 
cant. 

Figure 10 shows the two-dimensional finite Fourier transform 
of the map in Figure 5. A significant amount of the power is in 
wavelengths greater than 3000 km. However, we will argue that 
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Figure 10: Two-dimensional finite Fourier transform of a 

36° x 36° cen ered square subset of the bins of data o^ed 
construct Figure 1. Nyquist wavelength is 660 ku r.id the D.C. 
component is at the center. Contour interval is 40 nT/cycle. 
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Figure 12: Static field model produced by the two-dimensional 

finite Fourier transform filter described in Appendix II. 
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the source of a large fraction of this power is believed to 
reside outside the core. In Appendix 2 we describe an alternate 
method for deriving a scalar magnetic anomaly map using two- 
dimensional Fourier filtering rather than one-dimensional polyno- 
mial fitting to e«:cn pass. Figure 11 is the scalar anomaly map 
derived by this pr cedure and Figure 12 is the model of fields 
derived in the cor » (which have not been filtered by MGST(4/81)) 
and in the mjgr eiosphere. Figure 12 is just the difference 
between a *• • of B-BMOD and a high-passed version of this map, 
where b * ► X + + 2 ’?■ . Since the core field is relatively 

static over the five-month data window, Figure 12 is believed to 
be a good model of the fields generated in the core greater than 
4000 km in wavelength remaining after the removal of MGST(4/81). 
Since the qualitative correlation between Figures 5 and 11 is 
high and core field effects due to wavelengths greater than 
4000 km are believed to be successfully filtered from Figure 11, 
the long-wavelength power in Figure 5 is probably not due to the 
fields arising in the core which have not been successfully 
filtered. This argument only holds for wavelengths greater than 
about 4000 km. However, Langel and Estes (1982) have shown that 
the maximum amplitude of the core field between about 3000 km 
(GSFC(9/80) degree 13) and about 4000 km (GSFC(9/80) degree 10) 
wavelengths is approximately 20 nT at the earth's surface. (The 
core field model in our study removed terms of these wavelengths. 



Figure 13: Map of perturbations in the scalar field due to the 

horizontal anomalies found in Figures 6 and 7, if these anomalies 
are stable over MAG SAT's lifetime. 
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Figure 14: Schematic drawing showing how variations in the 

horizontal field can create perturbations in the measured scalar 
field. XMOD, YMOD, and ZMOD are the components of the core-field 
model; HMOD is the amplitude of the horizontal model field; BMOD 
is the amplitude of the scalar core-field model; AX and AY are 
anomalies in the x and y directions, respectively; H is the 
horizontal component of the measured field and BH is the measured 
scalar field resulting from the core field and horizontal 
anomalies. Then B1-BM0D is the perturbation in the scalar field 
caused by the horizontal anomalies. 
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Errors in the model between wavelengths of 3000 and 4000 km, a 
fraction of the field at these wavelengths, should be insignifi- 
cant in relation to the fields generated in the crust. There- 
fore, it appears that the long-wavelength magnetic anomaly fea- 
tures in Figure 5 are on the average not due to fields generated 
in the core. 

Second, since the major source of magnet ospheric fields mea- 
sured near the earth is due to equatorial ring currents found at 
an elevation exceeding three earth radii, the d.c. component of 
the ring current would generate a very long-wavelength magnetic 
field over Antarctica. Furthermore, since high power terms of 
the magnetospheric field change only very slowly in relation to 
the 25 minutes it takes the satellite to complete a pass between 
55 °S and 55 °S over the pole (Report of the Kakioka Magnetic Observa- 
tory, 1981), magnetospheric fields generated by ring-currents 
should be well-modeled by a quadratic polynomial fit to the data 
over Antarctica. (Actually, in hindsight, two quadratics separate- 
ly fit to half of each pass would have created a better magneto- 
spheric model). Magnetospheric fields should not find their way 
into Figure 5 in any significant proportion. 

This leaves us with the question: do fields generated by 
field-aligned currents seriously affect the scalar anomaly map? 

This is a difficult question. 
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Figure 16: Average vector core-field model, YMOD, in the in- 

direction (tangential). Positive clockwise. Contours in 10^ nT. 
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Figure 17: Average horizontal component of the core-field model, 

HMOD. Contours in -10^ nT. 
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Figure 18: Avera^ vertical component of the core-field model, 
ZMOD. Positive ci_t of page. Contours in 10^ nT. 
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The vector core field model (XMOD, YMOD, ZMOD) in Figures 
15, 16, and 18 varies over Antarctica so that the relationship 
between the horizontal field and the vertical field changes dur- 
ing a flight pass. In areas away from the geomagnetic south pole 
the horizontal component of the core model field, HMOD in Figure 
17, will be a large fraction of the vertical component, ZMOD, 
making variations in the horizontal anomaly /h * (AX^ +• AY^), 
more likely to contribute to AB. However, auroral phenomena away 
from the geomagnetic pole are less intense so that it is rather 
unclear which areas are expected to be most affected by the field- 
aligned currents. Nevertheless, in all auroral regions, varia- 
tions in AZ due to field-aligned currents should contribute to 
perturbations in AB. However, as shown in Section II, high- 
frequency field-aligned current effects do not find their way 
into AB, though longer-period variations in AZ and AH may still 
contaminate the scalar map. If long-period field-aligned current 
effects exist which perturb all three components of the vector 
magnetometer, then the scalar anomaly would show field-aligned 
current effects also. The horizontal anomalies are the places to 
look for the presence of field-aligned current effects. 

If we assume that Figures 6 and 7 represent the magnetic 
fields in the North and East directions generated by the long- 
period components (> 5 months) of the time-varying field-aligned 
currents, then Figure 13 is the map showing the perturbation 
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these fie.ds have on the scalar anomaly map. The values on the 
map were calculated in the way illustrated in Figure 14, in which 
BMOD is the value of the scalar core field model and BH is the 
scalar field containing perturbations due to AX and AY. The 
values of BH-BMOD for each bin is what is found in Figure 13. 

The poor correlation between Figures 13 and 5 probably implies 
that the AX and AY maps are not good representations of the long- 
period field-aligned current effects. However, if we retreat a 
little and only assume that Figures 6 and 7 indicate the general 
location of the long-period anomalies if, indeed, tl. xist, 
then Figure 13 would be interpreted as flagging :as ,£ poten- 
tial perturbation in the scalar anomaly map. The notable conti- 
nental areas sensitive to perturbations in the scalar field then 
would be the Pensacola Mountains region, the region inland of 
Thurston Island in Marie Byrd Land, southern Enderby Land, and 
the Amery Ice Shelf region. These features all lie within the 
region between 10° and 25° of the geomagnetic south pole in the 
region of highest field-aligned current activity. However, too 
much should not be read into this since most of Antarctica lies 
between these latitudes. To the extent that we believe that long- 
period field aligned current effects exist and that Figures 6 and 
7 capture the general location of these effects, we have to be 
concerned about the possibility of fields generated by field- 
aligned currents intruding the scalar anomalies in Figure 1. 
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Further tests concerning the existence of fields generated by 
long-period field-aligned currents are continuing to be run. 

Even though sensitive areas exist in Figure 13, two factors 
contribute to the belief that the anomaly featured in Figure 5 
are mostly not the result of field aligned currents. First, a 
necessary characteristic of fields generated in the earth's crust 
is temporal stability over periods during which the inducing 
field and the thermal conditions of the studied area change 
little. Comparison of Figure 5 with the POGO map generated by 
Regan et al (1975), shows the good qualitative correlation be- 
tween the MAGSAT and the POGO scalar anomaly fields. The differ- 
ences are mainly in amplitude, due to the lower orbit of MAGSAT, 
and in short-spatial wavelength features, perhaps Hue to MAGSAT' s 
greater resolution, data reduction differences between the method 
of Regan et al and ours, or seme field-aligned current effects 
not detectable using POGO's scalar magnetometers. In general, 
the correlation is quite good, indicating that over the 10 year 
period spanning the POGO missions and MAGSAT' s lifetime, the fea- 
tures represented in Figure 5 are stable. This lends credence to 
the proposition that the source of the anomalies in Fig. -.ire 5 lies 
in the crust. However, it is still unclear if there exists a 
field-aligned current component stable over this time period, sc 
this argument cannot be taken too far. 


Second, however, the correlation between the magnetic anoma- 
lies in Figure 5 and known geologic features, especially in ocean 
ic regions, is remarkable. Thu is interpreted to mean chat the 
crustal field is significantly represented in Figure 5 even where 
the horizontal field is believed to perturb the scalar field. 
These correlations will be dealt with in detail in Section IV. 
Thus, the temporal stability of the anomaly features in Figure 5 
together with the remarkable correlation between the anomalies 
and known geologic structures, lead us to believe that field- 
aligned currents probably do not severely contaminate the scalar 
anomaly map except in a few isolated regions. 

In summary then, the sources of the anomaly features in Fig- 
ure 5 appear to reside mainly in the earth's crust since sources 
below and above the crust appear to have significant effect only 
in exceptional regions. 

B. Vertical Anomalies 

The vertical anomaly map found in Figure 8 is in remarkable 
agreement with the scalar m^p. This is as we would expect in 
high latitude regions in which the earth's magnetic field is near 
ly normal to its surface. The rapid variations in AZ observed in 
Figures lb and 4b appear mostly to be filtered in the averaging 
and data-discarding process. However, the amplitudes of the 
anomalies in Figure 8 generally are greater than those in Fig- 
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ure 5. It is unclear which map better represents the crustal 
field, but they are similar enough to make this matter relatively 
unimportant. Vertical anomalies in Antarctica are of approximate- 
ly the same amplitude as those found elsewhere in the world. 

As shown above, the horizontal anomaly maps in Figures 6 and 
7 result mostly from fields generated external to the Earth. We 
suspect that the vast majority of the power has field-alig’.ied 
currents as its source. Thus, the amplitude of the horizontal 
anomalies, found in Figures 6 and 7, can be generally considered 
to be a good approximation to the amplitude of the field-aligned 
current effect. Therefore, horizontal anomalies over Antarctica 
have much larger amplitudes than those over other regions. How- 
ever, it is unclear if thi3 is related to any real long-period 
field generated by field-aligned current;!, which would perturb 
the a B map, or just to uncompensated random short period varia- 
tions which are believed not to effect the scalar anomalies. 

Crustal ly-genersted magnetic fields over most of the studied 
region will have a small horizontal component due to the nearly 
vertical nature of the core field. However, in the low geomag- 
netic Latitudes north of the Antarctic Peninsula, the earth's 
field has a strong enough horizontal component to register a hori- 
zontal anomaly of crustal origin. Unfortunately, however, the 
contour intervals in Figures 6 and 7 are insufficiently refined 
to detect a crustal signal in this region. Continental remanent 
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magnet iza cion may affect the horizontal anomalies if the earth's 
field was nearly horizontal during the period of batholith solidi- 
fication. The map resolution problem together with the latitudi- 
nal stability of Antarctica during the last 400 million years 
(implying a vertical geomagnetic field) imply that a remanent 
effect probably is not noticeable in Figures 6 and 7. 
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IV. GEOLOGICAL INTERPRETATION 

The interpretat ion of MAGSAT data, though still in its in- 
fancy, should hold great interest to the Antarctic geoscientist, 
for here is the first coherent continent-wide data set with infor- 
mation about the Antarctic crust. Furthermore, MAGSAT magnetic 
anomalies appear to be highly correlated with known Antarctic 
geologic and tectonic features, especially in oceanic regions 
where the geology is simplest and best known. 

Oceanic magnetic anomalies are almost invariably associated 
both with basins (negative)(cf . Frey, 1982a) and spreading ridges 
(positive). (All oceanic feature names will be taken from Heezen 
and Tharp ,. 1980 . ) Three of the four major oceanic basins surround- 
ing Antarctica (the Weddell, Enderby, and Wilkes Abyssal Plains) 
have strongly negative anomalies associated with them, and even 
the fourth (the Bellingshausen Abyssal Plain) is relatively low. 
(Geographic names are indexed in Figure 5.) The most striking 
conjunction of a positive anomaly with a spreading ridge occurs 
where the Mid-Indian Ocean Ridge and the East Pacific Ridge meet 
grid south of the Ross Sea embayment ("grid" directions refer to 
a Cartesian coordinate system laid across the polar map, with 
grid north parallel to the 0° meridian, grid east parallel to 
90°E, etc.). The set of positive anomalies running between 140 a E 
and 120°W, north of 65°S, lies closely over the East Pacific 
Ridge on the grid west and he Mid-Indian Ocean Ridge on the grid 
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east. There are also highs associated with aseismic volcanic 
ridges and plateaus, such as the Kerguelen Plateau (about 80 e E) 
and Maud Rise (65 # S, 0°E)(cf. Frey, 1982a), and a relative high 
in an otherwise pronounced low is associated with the South Sand- 
wich Islands and Trench (60°S, 25°W). 

There are, however, some interesting anomalies that do not 
fit this norm. For example, a positive anomaly runs grid north- 
east from Maud Rise right into the Enderby Abyssal Plain, and 
another extends grid west of Thurston Island into the Bellings- 
hausen Abyssal Plain. The cause of these anomalies is puzzling, 
and deserves further study. 

The correlation between magnet ic anomalies and continental 
structures is also striking. In East Antarctica, the mountains 
of Queen Maud Land (negative), the mountains of Enderby Land 
(positive), much of Wilkes Land (positive), the Gamburtsev Sub- 
glacial Mountains (negative), and the Amery Ice Shelf (negative) 
all have magnetic anomalies associated with them. Although it is 
not certain, of course, what these associations mean, some specu- 
lation may nevertheless be useful. We believe that the Enderby 
Land high may stem from a relatively high crustal magnetization — 
aeromagnetic surveys in parts of the area (Wellman and Tingey, 
1982) suggest to us that the mean susceptibility of the upper 
crustal rocks, at least, is higher than the continental norm. We 
suggest that the low over the Gamburtsev Subglacial Mountains 
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results from an elevated Curie isotherm; this idea finds some 
support in the low surface-wave group velocities along paths 
traversing these mountains (Dewart and TBksUz, 1965; see 
discussion in Bentley, in press), since a warming of the mantle 
causes seismic wave velocities to diminish. Perhaps these moun- 
tains are relatively young. The pronounced relative magnetic low 
overlying me .jnery Ice Shelf/Lambert Glacier region supports the 
belief that ..nis is a failed rift deeply filled by non-magnetic 
sedimentary rocks (Masolov et al, 1981); most continental rift 
features, unless characterized by extensive extrusive volcanism, 
show a negative anomaly (Frey, 1982a). The apparent extension of 
the anomaly into the ocean is probably just a failure to resolve 
closely adjacent continental and oceanic lows. 

On a larger scale, the coincidence of depressed topography, 
satellite-measured free-air gravity lows, and other features ex- 
tending from Wilkes Land across the ocean into Australia, led 
Veevers (1982) to suggest that the whole vast region is being 
held down dynamically by downward currents in the mantle. The 
positive magnetic anomaly in Wilkes Land (and the corresponding 
one in Australia — see below) is consistent with this suggestion 
since a convection-convergence zone would be relatively cool. 

A noteworthy feature of the anomaly map is the absence of 
magnetic anomalies over the Transantarct ic Mountains. Instead of 
exhibiting a characteristic anomaly pattern of their own, they 
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mark a distinct boundary zone between largely separate East and 
West Antarctic anomalies. (It is likely that the negative anoma- 
ly that cross-cuts them from the Ross Ice Shelf is, in fact, 
another case of two separate lows that have not quite been 
resolved. ) 

Several anomalies appear in West Antarctica, but their 
tectonic association is not at all clear. The volcanic province 
or provinces comprising the Antarctic Peninsula, Thurston Island 
and Marie Byrd Land all show distinct highs as would be expected, 
but along the whole region the centers of the highs are inexplic- 
ably shifted oceanward. A pronounced low over the Ross Sea may 
support tne concept of a failed rift zone here, but it is not 
centered over the postulated axis of the rift found from gravity 
measurements in the grid eastern part of the sea (Hayes and 
Davey, 1975; Bentley, in press, Figure 9). The anomaly does 
disappear under the Ross Ice Shelf — that is in agreement with 
the gravity evidence (Davey, 1981). 

What may be a mirroring negative anomaly appears in the 
Weddell Sea embayment between (and partly overlying) the Ells- 
worth and Pensacola Mountains. However, if this is a rift-zone 
negative, it is surprising that it does not extend farther grid 
northward under the Weddell Sea continental shelf. 


Comparison of the Gondwana reconstruction of Norton and 
Sclater (1979) (Figure 19) with Figure 5 and the global anomaly 
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map of Langel et al (1982a) shows that the pronounced high in 
Wilkes Land is mirrored by an even more pronounced high in the 
Australian shield; there is, in fact, a general similarity be- 
tween the magnetic appearances of Wilkes Land and Australia. 
Moreover, the low in Queen Maud Land appears to correspond well 
to lows in southern and southeastern Africa. On the other hand, 
there is no clear correspondence between the West Antarctic posi- 
tive anomalies and any g else. It appears, unfortunately, 
that the MAGSAT map is not yet going to solve the puzzle of where 
to put the West Antarctic microplates before the break-up of 
Gondwanaland! For East Antarctica, nevertheless, the Gondwana 
magnetic reconstiuction is very good, implying that the anomaly 
features in Antarctica were formed prior to break-up. A more 
exact comparison requires that all data be reduced to the pole. 


-56- 


V. CONCLUSION 

We have argued that of the four high-power sources of mag- 
netic fields in high latitudes (the core, magnetospheric ring- 
currents, field-aligned currents, and the crust) the crust is the 
source of most of the anomalies observed in the scalar (Figure 5) 
and vertical (Figure 8) vector anomaly maps and that field- 
aligned currents generate the anomalies seen on the two horizon- 
tal anomaly maps (Figures 6 and 7). 

Errors in the core-field model do not provide a significant 
fraction of the power observed in Figure 5. Errors in the model 
contained in degrees 1 and 2 are compensated by the quadratic 
polynomial fit to each pass, and errors contained in degrees be- 
tween 10 and 13 are known to be too low in power to seriously 
contaminate Figure 5. Errors in terms between degrees 3 and 9 
may appear on the scalar anomaly map, but the power between these 
wavelengths is not a large fraction of the total power of the map. 
Fields generated by magnetospheric ring-currents are known to 
vary slowly relative to the time it takes MAGSAT to traverse the 
observed region and are therefore well modeled by the quadratic 
fit to each pass. Finally, though the horizontal anomaly data 
indicate that field-aligned current activity is significant on 
every pass and that high-frequency variations in these currents 
cause the anomalies seen in the horizontal anomaly maps, we do 
not believe they are significantly contaminating the scalar anoraa- 
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ly map. Short-period field-aligned current effects do not dis- 
turb the scalar anomalies, though longer-period currents can. 
However, very long-period field-aligned current effects have not 
been observed, not because they do not exist but because they 
must be much lower in power than the high-frequency variations. 
The good correlation between the scalar anomaly map and the POGO 
anomaly map over Antarctica implies that the source of most of 
the power of the map is stable over 15 years. Since field- 
aligned current terms longer than 15 years in period are not 
strong enough to cause the observed field, long-period field- 
aligned currents are believed not to significantly alter the 
scalar anomaly map. Furthermore, the remarkable correlation of 
magnetic anomalies with known oceanic geology implies a crustal 
origin of the oceanic anomalies even in ' r ery high geomagnetic 
latitudes . 

The scalar anomaly map is consistent with the following hy- 
potheses that we associate with the cited authors: 

1. The Gamburtsev Mountains are tectonically unrelated to 
either the Highland Massifs (Drewry, 1975) of Wilkes 
Land or the mountains of Enderhy Land; 

2. Wilkes Land is the site of cool convergence in the 
mantle which accounts for the topographic depression and 
gravity low (Veevers, 1982); 
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3. West Antarctica comprises the set of tectonic micro- 
plates postulated by Dalziel and Elliot (1982); and 

4. Both the Ross Embayment (Hayes and Davey, 1975) and the 
Amery lee Shelf mark the site of ancient continental 
failed rifts. 

Finally, a comparison of the anomalies observed over Antarc- 
tica with anomalies over the other Gondwanaland continents indi- 
cates that the source of the anomalies on the scalar map has been 
stable since before rifting of Gondwanaland occurred. 

In the future, we will attempt a comparison between existing 
aeromagnetic surveys in Antarctica and the MAG3AT results. The 
only data available to us are from West Antarctica (extensive 
work has also been done east of the Weddell Sea by the Soviet 
Union, and in the Enderby Land region by Australia). It remains 
to be seen whether the lateral extent of the West Antarctic sur- 
vey is great enough to yield a meaningful upward continuation to 
MAGSAT heights. 

Modeling of MAGSAT anomalies is of questionable value at 
this stage, for two reasons. First, too little is known about 
the large-scale structure of the anomaly-associated features in 
Antarctica to provide any realistic geologic control on the 
models. Second, we believe further resting for auroral-zone 
current effects by comparing dusk-side passes with dawn-side 
passe; is needed. If it can t’ en be shown that the data are 
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sufficient ly decontaminated, we would suggest looking first at 
the Queen Maud Land anomaly, together with the associated nega- 
tive in Africa where the geology is better known. The anomalies 
in the Enderby Land/Amery Ice Shelf area and in Wilkes Land will 
presumably be attended to by MAGSAT principal investigators from 
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Abstract . A procedu. iclecc and rsducs 

•ateilics u^itcic anoea. uses La high southern 
latltudts Is dsscrlbsd, and a map of Antarctica, 
constructed using this procedure, is snovn. The 
nap is qualitatively analyzed for error and geo- 
logic significance. Correlations are noted 
between magnetic anonelles and aountain ranges, 
subgiacial basins, tectonic provinces, regional 
gravity anomalies, a hypothetical continental 
rift feature, oceanic basins, and oceanic rises. 
Overall, the correlation between the negnetlc 
anomaly patterns and known geological features 
la good. 

Introduction 

The Earth's magnetic field is principally a 
superposition of fields generated in four 
regions: the core, the crust aoove the Curie 
isotherm, :ne ionosphere, and the magnet ospnez «. 
Isolating tee crustal field at satellite 
elevations is difficult, especially at high 
latitudes, though attempts span sore than a 
decide Lc.g. Celt:, et el.. 1970; Langel, 1974; 
Regcn. et al. , 1975; Coles, et al. , 1976; Coles, 
1979; Coles. «c al. , *979; Mayhev, et al., 1930; 
Langel, et al. , 2980a). The difficulty arises 
due to the proximity to Che satellite of the 
ionospneric and magnetospnerlc currents, vnich 
are amplified in auroral regions. In this paper 
we describe a aechod for extracting the crustal 
field frea the total field aeasured oy MACS A I 
at high latitudes and give a brief interpretation 
of some features of the Antarctic crustal eag- 
nctic jnctaaly field. 

Data Selection and Reducclco 

The crustal aagnetic anomaly, <18, Is cal- 
culated by subtracting a core field aodel and 
an external field aodel freo the total scalar 
aagnecic field aeasured by MAGSAT. The degree 
and order 13 spherical naraonlc nodel MCST-/8L 
created by Langel. et al. (1930b) wee used as 
the cere field zodel. No such general external 
field mcael exists; therefore, ve followed the 
procedure described below in filtering external 
field effects. 

External fields primarily consist of two com- 
ponents: a long-wavelength, slow-varying field 
generated by ring currents in the aagnecosohere 
and a snorter wavelength, faster-varying field 
generated by fie Id- aliened currents in cae iono- 
sphere. The rapid variations in field-aligned 
currents soke then difficult to aodel. so ve 
deleted all passes snowing signs of their effect. 
This was done in two ways. First, we selected 
data only from passes during which the planacary 
magnetic activity index, Cp, was leas than or 
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equei to 1“ for at least 6 hours. The extended 
tlae period wea chosen because auroral lor.o- 
spnerlc activity abatei acre slowly than low- 
lacltuds activity. Second, data wera saiacted 
only from passes showing no sign of anomaly 
amplitudes greater than thoaa theoretically ex- 
pected from the crust. This was done since 
fleld-ellgned currents, because of their 
proximity to tho satellite, can generate much 
larger aeasured magnetic effects than tha cruse. 
Specifically, wa have required chat IB be no 
larger than 15 games (nT). This amplitude 
criterion can be applied only after the long- 
wavelength external field is modeled and removed. 
Over the 8000 Icn flight path ve have considered, 
the long-vevcicngth field ij modeled as a Quad- 
ratic polynooial (a.g. May new. 1979) upon which 
the crustal and ionospheric fields are super- 
imposed. 

An example of this procedure is given in 
Figure 1, in which traces A result froa the 
meesured magnetic field after the core field has 
been removed, traces B ars tbs models of the 
field generated by the ring currents, and traces 
£ depict -3, which results by subtracting 3 from 
A. Tha trackj of tha 37 passes (cut of a total 
of approximately 2300) over Antarctica between 
November 1, 1979 and April 1, 1930. that 
satisfied the selection criteria are snovn in 
Figure 2. 

Tha IB data from the selected passes were 
averaged over arees measuring 3* of latitude by 
3* of longitude, and were plotted and contoured, 
yielding the scalar anomaly aao snovn in Figure 
3. Anomalies cn this map have not been 
corrscted for altvatlon variations in tha 
satellite path, nor have they been reduced to 
thr pcLe. 3och correctiooe will eventually be 
applied, since the 200 km elevation variations 
that occur between MAGSAT orbits alter tha 
amplitudes even of very long wavelength aagnetic 
anomalies by as much as a factor of two (Regan. 
1979; Bhattacharyya, 1977), and since the geo- 
magnetic latitude varies by more than 45* across 
tha map so that fields induced la the crusc by 
the Earth's dipole field are lacltudlnally 
dependant. At presenc, however, we assume char 
elevation variations evenly sneer aoec anomalies 
and not# that for most of tha map, reduction to 
tha pole will have amimal effect. Ocher re- 
finements of tha data planned for tha near future 
include rectangular grlddlng, tha incorporation 
of a cross-correlation selection criterion for 
adjacent and crossing satellite paths, and data 
selection based on analysing the vector -magnetic 
rnomaly data for fie id- aligned current signatures. 

Accuracy of Data 

Disregarding magnetometer and tracking error, 
the accuracy of the magnetic anomaly map La 
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Fig. 1. Craphs of aagneclc anomalies for: 

(a) a past unaffected by field-aligned currents, 
and (b) an arxecced pass wrier « A is the total- 
field ainus cne core field aidel, 3 is the quad- 
ratic polynomial aodel of the aagnetospncrlc 
field, sod C is A ainus 3. i.e. IB. 

Figure 3 depends on the degree to which the core 
field aodel and the external field aodel are in 
error. We concerned ourselves chiefly with 
errors found in the external field aodcls, since 
they are both aore primitive and shorter in 
wavelength than the core field aoael. Two tests 
based on the tiae-varying nature of the external 
fields were performed. First, we partitioned 
the MACS AT data into two successive Z^-aontn 
data secs and eaployed the data selection and 
reduction schese described above to each. The 
correlation between the two aaps was good, 
suggesting that most of the external field 
effects with periods less than five months were 
successfully filtered. Second, we c smeared the 
MACSAT map with the PCCO nap of Regan, cc al. 
(1975) covering Antarctica. The same general 
continental features exist on both maos, further 
suggesting chat sost external field effects with 
periods less than ten years were filtered out. 

The results of these tests are encouraging, but 
a aore telling test, previously oerforaed for 
high northern latitudes oy Langel. et al. (1980a) 
and Coles (1979), would be to compare Figure 3 
with lover altitude data continued up to 
satellite elevations. Until rcccncly. data were 
not available froa any aagnetlc survey over 
Antarctica that waa large enough la area and 
dense enough in coverage to warrant comparison 
with MACSAT. However, an aeromagnetlc aao soon 
to be published (referred to bv Janxovskl. et 
al. , 1981), covering alaoec 10° <■- in West 
Antarctica, will provide an excellent check of 
MACSAT lata. 


Source of Anonalles 

7e expect that two aajor factors will con- 
tribute to continental cruatal aagneclc 
anomalies: Che depth to Curie isotherm end the 
magnetic susceptibilities. It is still not 
certain to what degree continental remanent 
magnetisation may affect long-wavelength mag- 
netic anonalles. Some workers, at lease, con- 
sider the efface negligible iM. Mavhev, personal 
communication, 1981); lacking better information 
we Ignore it here. Thus we assume chat the main 
sources of aagneclc anomalies are located above 
tha Curie Isotherm and above Mono (Wasilevskl, 

•t al., 1979), and principally in cne lower 
crustal layer where susceptibilities are greatest 
(Hall, 1974). Therefors. continental hlgha in- 
dicate either a thick crust with a deep Curie 
isotherm (low heat flow), or exceptionally high 
aagneclc susceptibility in chc lower crust, or 
both. On the other hand, continental Lowe imply 
some comoinat ion of a thin crust, a shallow 
Curie isotherm (high heat flow), and low 
susceptibilities 13 the lower crust. 

The highest susceptibility region in the 
oceanic crust is id Layers 2 and 3, nearer the 
surface chan in cne continental crust, and 
generally veil above the Curie isotherm. There- 
fore, oceanic magnetic anomalies should reflect 
regional crustal thicknesses and susceptibility 
differences, but probaoly do sot reflect heat 
flow except in extreme circumstances. Remanent 
magnetization In the rocks of oceanic rift zones 
can also yield magnetic anomalies at satellite 
elevations, as can be snovn oy model calcula- 
tions. 



Fig. 2. Flight eracks of the 37 \ccepced passes 
over Antarctica and the surrounding oceans. 
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Geological Interpretation 

The Interpretation of aatcllica magnetic 
Anonaiy data la clterly ambiguous, especially 
ov«r regions with Uetla collaearal data, iov- 
•ver, tha major associations between geologic 
features and magnetic anomalies ara worth noting 
la both cootloaotai and ocaanic raglona. 

Ovar tha Antarctic conciaanc chart ara 
■agnatic anomalies aaaociatad with zany topo- 
graphic and gtologlc faacuras. In last 
Antarctica, tha mountain* of waatarn Queen Maud 
Land (magnetic Low), tha aountaina of Endarby 
Land (high), tha ?rlnca Char las Mouncalna (low), 
rha "Highland Massifs" baaaach tha ica In Wilkes 
Land (high) (Trevry , 1975), tha Gamburtsev Sub- 
glacial Mountains (low), and tha Pensacola 
Mountains (high) all have magnetic anomalies 
aaaociatad with them. Also, thara is a magnetic 
low ovar a subglacial basin generally along 150*1 
in East Antarctica thac is furtnar charactarlzad 
by a strongly negative regional lsoscatlc gravity 
anomaly. This suggests that tha process associ- 
ated with the lsoscatlc imbalance say ba associ- 
ated with hign neat flow and does not support a 
nodal cf cool convergence in the circulation of 
tha upper mantle. 

In West Antarctica, the Antarctic Peninsula, 
eastern Marie Byrd Land, and the Thurston Island 
region ara associated with magnetic highs, where- 



Fig. 3. MAGSAT total-field magnetic o.nomaly sap 
ovar Antarctica, '.’nits in gasses mT). Capital 
letters indicate the approximate location of: 

A - Quean Maud Land. 3 - Enaerby Land. 

C - Prince Charles Mountains, 3 - Anar lean 
Highlands, E - Canourtsev Mountains. F - Wilkes 
aid, G - Transaotarccic Mountains, H - Rosa 
ibayaaac. I - Maria 3yrd Land. J - Thurston 
■land, K - Ellsworth Mountains, L - Pensacola 
Mountains , M - Antarctic Peninsula. 


aa the Ronna tea Shelf /Ellsworth Mountains region 
shown a pronounced aagnatic low. Several of 
chesa anomalies seta to be closely Associated 
with tacconlc alcroplacas suggested by 3.H. 

Ellloc (personal cconunlcatica. 1981). Hayes 
and Davey (1975) have argued from tha occurrence 
of a striking linear positive gravity anomaly 
that an actual ancient rift cone, whlcn was 
probably active during tha separation of 
Antarctica from Australia 65 m.y. ago, underlies 
che western Rosa S as. If so, the crust snould 
be anomalously thin and tha hast flow higher 
than could be normal there. * magnetic low 
covers such of this area, although the trende of 
the two feacurca iiffer. 

Tvc continental regions thac differ only in 
crustal thickness will tand to show different 
magnetic anomalies, with sore positive values 
over tha thicker crust. 3enclay (1973) iod 
Dcwarc and Toxsoz (1954) have estimated crustal 
chicanes* to be 40 km in East Antarctica ana 
30 tan in West Antarctica; tha nagaecic anomalies 
over East Antarctica are generally sore positive. 

Comparing the Conavana reconstruct ion of 
Craddock, ec ax. (1970) with Figure 3 and tha 
global magnetic anomaly sap of Langel, ec al. 
(this issue), shows thac the Low over tha Ross 
Sea and the Transanearctlc Mountains corresponds 
to tha low along the Adelaide and Tasnan orogens 
in eastern Australia, that tha hign in Wilkes 
Land is mirrored by a high in the Australian 
shield, and chat, in fact, cnere Is a general 
similarity between Che aagnetlc aopearances of 
Antarctica and Australia. Moreover, the Enderby 
Land high seems to be expressed in che Indian 
Shield, though this is unclear cue to the 
latitude difference between tha two region*. 

Tha Queen Maud Land low appears to correspond to 
a low m southern Madagascar *nd cne Ellsworth 
low mry be reflected in a low over tha Case 
orogen in South Africa. On tha ocher hand, co 
feature similar to the Antarctic Peninsula hign 
is apparent along the Andean orogen in South 
.America. Thus, the Gondwana reconstruction 
seems largely consistent Vith the magnetic data 
snown In Figure 3 and Langel' s mao, but some 
discrepancies exist. A more exact comparison 
requires reduction of all data to the poLe. 

Oceanic magnetic anomalies ire associated 
with both basins and rises. All four major 
ocean basins around Antarctica exhibit magnetic 
lows, and with the exception of the linear mag- 
netic low along 170*V, all the oceanic negative 
anomalies are associated with these basins. 
Oceanic basin anomalies presumably are mainly 
negative because of the tnln oceanic crust. 
Conversely, many of th. oceanic rises in the 
region exhibit higns, and note oceanic magnetic 
highs 11a above known rises. However, not all 
oceanic rises exnlbit a magnetic high. We 
relieve che prinarv cause tor the magnetic hlgna 
Is remanent ragntcizacion over spreading ridges, 
althougn soeciflc modeling has noc yet been 
carried ouc. 

Conclusions 

Preliminary tests indicate that che primary 
source of che magnetic anomalies snown i *igur* 

3 lies within the Earth's cruse. Thougn tha 
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data used to cooatruct the magnetic anomaly uo 
were aw: corrected for the effects of elevation 
*od latitude varlatioaa in the satellite flight 
path, both contlneatal and oceanic anomaly f»a- 
turea ihov a good agreement with known geologic 
structures. Further research will focus on im- 
proving the quality of the magnetic anomaly map 
and metaling the sources of the anomalies seen. 
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At -NDIX II 

Alternate Data Reduction Procedure: 

Two-Dir >nsional Finite Fourier Transform Method 

Filtering polynomials from MAGSAT passes in non-polar 
regions where passes do not cross is the source of the latitudi- 
nal striping apparent in Langel et al (1982). Removal of the 
polynomials amounts to removal of most long-wavelength North- 
South power so that only long-wavelength East-West power remains; 
thus, the East-West striping. In polar regions, flight paths 
cross so that this problem is not generated. Therefore, Antarc- 
tica is the perfect proving ground for testing alternate data 
reduction procedures to the standard procedure employing polyno- 
mial fitting. Cain et al (1982) have tested a whole-earth two- 
dimensional spherical harmonic method with what they feel to be 
surprisingly good results. We have tested another method on a 
continental scale with somewhat ambiguous, though heartening, 
results. 

Instead of subtracting a polynomial from each pass we have 
merely taken the uireduced data (observed values minus core-field 
model for each pass) and averaged them in 3° square bins. The 
resulting map was then high pass filtered using a two-dimensional 
finite Fourier transform (2DFFT) filter so that spectral peaks 
between 4200 km and 5280 km were diminished by 1/3, peaks corre- 
sponding to wavelengths greater than or equal to 5280 km were 
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dirainishe^ by 1/2, and the d.c. component was set equal to zero. 
The surface subtracted in this way is shown in Figure 12. This 
surface differs from the polynomials used in the standard method 
in that it is static over the 5 month data window whereas the 
polynomials in the standard reduction procedure can be viewed as 
two-dimensional slices of a three-dimensional surface varying 
during the same 5 months. The map produced can be seen in 
Figure 11. The question is, how well does the static model 
correspond to the dynamic model? 

The difficulty in comparing the effects of the two models is 
exacerbated by the fact that the no-data region due to the tilt 
in the satellite's orbit had to be filled in by linear interpola- 
tion to perform the Fourier transform. The spectrum will be 
affected, probably by adding long-wavelength power, since the 
linear interpolation acts to create a long-wavelength feature. 
Furthermore, the zero-levels of the two maps are not the same — 
they are the averages of each map, separately, and since Figure 
11 does not possess as much oceanic region as Figure 5 its aver- 
age will be somewhat higher — approximately .4 nT. Thus, fea- 
tures on Figure 11 are d.c. shifted by about .4 nT below features 
on Figure 5. Comparing the maps indicates that there is a pretty 
good qualitative correlation — highs in Figure 5 are matched 
with highs in Figure 11, and lows with lows. There are two major 
exceptions, the Antarctic Peninsula high in Figure 5 is matched 
to a low in Figure 11, though this region is still a high reia- 
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tive to Ellsworth Land and the surrounding oceans. Also, the 
radial feature found at approximately 75°S and 150°W is not found 
in ''igure 5. In general, however, the qualitative agreement be- 
tween the two maps is good. The radial striping in Figure 11 is 
reminiscent of that found in Ritzwoller and Bentley (1982), which 
was believed largely to result from field-aligned current effects 
disturbing individual passes. 

The quantitative agreement is not so good, however. Ampli- 
tudes can be a factor of two greater in Figure 11 than in Figure 

5. This is the result of the fact that the static model cannot v 

accomodate daily and seasonal changes taking place in the magneto- 
sphere or long-spatial wavelength changes due to field-aligned 
currents. The latter is a problem specific to high latitudes. 

Keeping the problems specific to polar regions in mind, 
i.e., anomalies in Figure 11 can be generated by field-aligned 

j 

currents and the data gap, the correlation between the two maps 
should be considered heartening to the low-latitude investigator 
who mourns for North-South r netic features. 

In conclusion, experimentation with modifications of the 
static field model by low-latitude investigators is believed to 
be warranted. Two-dimensional filtering should be carried out 
after the removal of a magnetospheric field model (e.g. ‘he model 
of Langel et al (1982a)). If the magnetospheric field model 
contains fields generated solely in the magnetosphere, subtrac- 
tion of this model will not create a striping problem. Perhaps 
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then magnetic anomalies associated with such North-South features 
as the Mid-Atlantic Ridge, the East Pacific Ridge, and the Andes 


may appear more clearly. 
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APPENDIX III 
Core Field Model 

Langel ... al (1980) describe Che procedure by which a geomag- 
netic core field model (MGST(6/80)) is constructed from the 
MAGSAT vector data. A similar procedure is followed in the con- 
struction of the reference field used in this study. The geomag- 
netic field is expanded in terms of spherical harmonics such that 
the potential, V, from which the geomagnetic field is derivable 
is given by: 

N n 

V = a l (a/r) n+ ^ V (g m cos m<t> + h m sin m<J>) P m (Q) 

— 1 A ^ ^ ^1 

n=l m=0 

where a denotes the radius of a reference sphere, r the geocen- 
tric radial distance, 0 the geocentric colatitude and <{) the longi- 
tude of the point at which the potential is to be found, and 

P m (0) d'-uote' • associated Legendre function of degree n and 
n 

order m. 

A complete description of the earth's geomagnetic field 
should extend to n = m = °°. However, Cain ( i9 76 ) has argued that 
a thirteenth degree and order representation is sufficient to 
represent the main field of the earth, and MGST(4/81) contains 
just these terms which ere similar to those listed by Langel et 
al (1980). The average vector core model field (XMOD, YMOD, 

ZMOD) is shown in Figures 15, 16, and 18. The figures were separ- 



ately constructed by taking the sector model values at each data 
location and elevation and averaging in the same 3° square bins 


used for the construction of the anomaly naps. Thus, 
figures represent the average vector core field model 
over the 5-month period of data acquisition. 


these 

employed 
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APPENDIX IV 

ExCernal Magnetic Fields 

This discussion is largely based on the excellent review of 
the current patterns external to the earth written by Schunk and 
Nagy (1980). The regions of magnetic field generation relevant 
to MAGSAT research lie in the magnetosphere and the ionosphere. 

The magnetosphere is that region where forces generated by the 
earth's magnetic field dominate and extends to approximately 
9 earth radii. The ionosphere is the layered region of near- 
earth plasma lying between about 90 and 500 km elevation formed 
in the suniit hemisphere by ionization of the earth's tenuous 
upper atmosphere. A schematic illustration of the earth's bow 
shock and magnetosphere is shown in Figure 20. 

Sunward, shocked solar wind plasma in the magnetosheath can 
flow directly into the ionosphere via the polar cusp. At MAGSAT 
elevations the cusp occupies a narrew latitudinal band centered 
near noon with considerable longitudinal extent. Within this 
band, energetic particle precipitation can provide measurable 
magnetic disturbances, but the disturbances am less appreciable 
than field-aligned current effects. 

The magnetospheri c n.antle is an intermediate region through 
which magnetosheath plasma travels on its way to the plasma sheet. 
The ring curren; represented in F Lg ire f.l if ; he earthward exten- 
sion of the plasma sheet. Trapped energetic protons and elec- 



Figure 20: Schematic illustration of the earth's bow shock and 

magnetosphere. Taken from Schunk and Nagy (1980). 



fig ure 21: Schematic representation of a partial ring current 
(bold arrows) and field-aligned currents (thin arrows). Taken 
from Anderson and Vondrak (1975). 
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crons grad-B drift in opposite directions causing ring currents 
to flow around the earth. The current flow is large enough to 
effect magnetic measurements greatly and is believed to be the 
major source of level bias variations between MAGSAT passes over 
Antarctica. 

The electric potential difference across the tail of the 
magnetosphere pointing from dawn to dusk, generated by the inter- 
action of the shocked solar wind with the geomagnetic field, is 
mapped into the high-latitude ionosphere as an electric field 
normal to the geomagnetic f^eld. This electric field generates 
the auroral electrojets, flows of ions parallel to the earth. 

Only the high-latitude ionosphere is influenced directly by the 
magnetospheric field since the ring-current effectively shields 
plasma from leaving the plasmasphere at low latitudes. 

The horizontal ionospheric field is further coupled to the 
magnetosphere through field-aligned (Birkeland) currents. The 
field-aligned current pattern during a quiet and a disturbed 
period is shown in Figure 22. The current patterns are concentrated 
in two areas encircling the geomagnetic poles, one poleward and 
one equatorward. The current flow is into the ionosphere in the 
morning sector and away from the ionosphere on the evening- sid»» 
in the poleward current region, whereas flow is opposite to this 
at a given local time in the equatorward regions. The net flow 
appears to be inward in the morning sector and outward on the 
evening side in the northern hemisphere. It is the field- 
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Figure 22: Summary of the d' ^ribution of flow directions of 

large-scale f ie ld-a lign<*'’ currents for (a) weakly disturbed condi- 
tions and (b) active periods. Taken from Potemra (1979). 
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aligned currents that are believed to be the principal cause of 
magnetic disturbances observed in the MAGSAT data. 


3 


APPENDIX V 


Computer Programs 


List of Programs: 


1. IBMTR 

2. REFMT 

3. ANTAP 

4. MRGRD 

5. BIN 

6. RDBIN 

7. FFT2D 


/ 
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O l ' 

Oh i (. v'.i '- v . i v 




C RlbMTK TRANSLATES THE I HM-RORMATT KIl NASA 1 APR 

C TO HAHHlS-RORflATTEH COM. THE KEY TO THE PRO- 

ORAM LIES IN THE SURHOUUNRS S. K'IRNSl. AMI S. 1TRNSL. 
C WHICH 1 RANSLATF. THE REAL AMI 1 HE INTEGER WORDS 

C RESPECT JVRI.Y. THE OLITPUT IS RORMATRJ) NEARL Y 

C IDENl 1CAU.Y T(J 1 HE E OH MAI GE 1 HR NASA 1 AI S E. 

C EACH HASS HKilNS A1 VOS W11H A HEADER RECORD. 

$ ROLL. OWE) i HY A STRING C)E DATA RECORDS CONTAINING 

C INE'ORMAl ICiN AHOLIT SO MAGNETIC VALUES. 1 HE 

C HEADER RECORD IS EORMATEH AS HOLLOWS: 


C 

C: WORD ft 

C 1 

C: 2 

C 3 

C: 4 

t: b-A 

C 7-8 

C V— 10 

f: 11-12 

C 13-14 

C 1S-1A 

C 17-22 

C. 23-34 

C 3b 

C 3A 

C 37 

C 38 

C 3V 

C 40-328 

C: 32Y-S24 

C b2ti-S27 

f: 

C 1 HE DA I A RECORDS 

C 

C WORD 0 

C 1 

C 2 

C 3 

C 4 

t: & 

C A 

(: '/-:«<* 

t: 37-AA 

C. A7-VA 

C: V7-1VA 

C 12/-JbA 

( lb'/- ISA 


NASA VamjAHLE 
I TYPE X 
NTYPEX 
MJUX 
I PASS X 
ASCX( 2 ) 

DSCXt 2 ) 

MSECXL 2 ) 

ALT MXC 7 ) 

AL.ONXt 2 ) 

1KP( 2 ) 

GSM( 7.i 3 ) 

DST 17. A ) 

NMAX 
NMAXT 
MODE XT 
T7ER0 
ASAR 

GH( 17. 17 ) 

GHT ( 14. 14 >2b 
fc{ 3 ) 

ARE FORMATED AS FC'LI OWS: 

NASA VAK1AKI.E 
I TYPES 
NTYPEH 
HODS 
MSECS 
1PASSR 
TINTS 
l.AT 
ION 
RAH 
Ml 1 

I NVI A I 
DIM A1 


(2 


one? jo ►- onoonnoooooooorjooooon 
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1H7-716 

HS 

217-246 

BV 

247-276 

X 

277-306 

Y 

307-336 

Z 

337-366 

HVA 

367- 3 V6 

XA 

W/— 426 

YA 

4I2/-4S6 

7A 

4M7-406 

BVSP 

4B7-S16 

XSCI 

S 17-546 

YSD 

54/-S76 

zsn 

S/7-606 

BMD 

606-636 

XMD 

637-666 

YMD 

66/-6V6 

7 MU 

6V7 

GOAL. 

6VH 

SPARE 


FOR A DESCRIPTION OF THF. NASA VARIABLES 
SEE. "HAGSA'I INVEST 1GAT0R-B TAPE GENER- 
ATION PROGRAM REQUIREMENTS*. NASA PUB- 
LICATION. CONTRACT NASS— 243V J . JULY, 1VGO. 


««•«•< «««*****♦•« «-**■*•«*#*»*»* ***■***■***■•#-**#*•**#**** 

INTEGF_R*6 NF.WVAR! 40 ) 

DOUBLE PRECISION NEWORD! 2 AO ). DX< 40 ). HARRIS! 7&0 ) 

INTEGER HUK 1010 L PC 3030 ). I COUNT . I WORD. RCOUNT 
CALL LINK 
WRITE! 1.1) 

FORMAT! ' FOR INPUT TAPE' ) 

CALL ASSTAP! ' 12 ) 

WRITE! 1.2) 

FORMAT! ' FOR OUTPUT TAPE' ) 

CALL ASSTAP! '13 ) 

•THIS LOOP CREATES AN ARRAY CALLED » UHICH IS 
COMf-’R) '/FT) OF WORDS EACH POSSESSING ONE BYTE OF* THE CORRES- 
PONDING WORD FROM THl AKRAY NAMED HUE! J ). THE BUFFER IN 
C: HEFOIsF. IT Ell. IS BUT. 

LOOP 

BUFFH.it IN! '12. HUE. B. 1010. 1ST Al . ’WORD ) 

CALL. STATUS! '12 ) 

EX) II OOP IF! 1ST AT. EG. 3 ) 

FOR J»J , )WOKD 
I 

B! 1 )=HUF!.I). AND. '77600000 
B! I )«H! ) ) SHIFT. -16 
M! !♦) )*"HUF! , I ). AND '177400 
FK 1 ♦ 1 )« M< 1 ♦) ). SHIFT. -H 


uor > c o 
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I Vi - . . y 

* k 

Cf V'- • - 


H( )*y )*WIK( J ). ANH. '377 
F.NPF OK 
C: 

(:«*««**TKNST..A I ION OK 1BH KF.AI..S INTO HARRIS KFALS 
C 

i can. ini - i 

RCOUNT M 

] F! IWOKR KC4. 1006) GOTO 6 
C 

(.•••«*»7KANSI.A7 ION OF HEADER RECORD 

CALI . 1TRNSL! 1 COUNT, 1. 16, B. NFWVAR ) 

CAl I RTRNSL! RCOUNT. 17. 32. B, NEWORD ) 

CALL 1TRNSI ( I COUNT . 33. 40. ». NFWVAR ) 

CALI. RTRNSL! K CXftJNI . 41 , 56. B. NGW0R1I ) 

CALI ITRNSL! ICOUNT. i>7. 64, B. NFWVAR ) 

CALI RTRNSL! RCOUNT. 65. 136. H. NI-.WORn > 

CAl 1 ITRNSL! ICOUNT. 257. 266. B. NF.WVAK ) 

CAl J.. R7RNSI ( RCOUNT. 269. 2228. B. NF.WORI) ) 

OCJTO 300 
C 

C««*-- «« TRANSLATION OF DATA RFCCRD 
A CAIJ_ JTRNSL! ICCUNT. X, 20. R. NEUVAK ) 

CALL fORNSL! RCOUNT, 21. 2784. B. NEWORD ) 

CALL ITRNSL! ICOUNT, 2785. 2904 . R. NF.WVAK ) 

CALL RTRNSL! RCOUNT . 2905, 3034. B. NKWORH ) 

•CREATE A NEW PBL PRECISION KtAI. AT<AY HARRIS FOR BUFFFROUT 

■•FIRST. UFLOAT NEWVAR 
FOK 1-1. 40 
BX( I 5-NFWVAR! I ) 

FNllFCJR 

IF< 1WGRD. ECL 1008 ) THEN 
FOR 1-1,5 
HARRIS! 1 >-UX! I ) 

FNIFOK 
FOK 1 -A. 696 
HARRIS! 1 )-NFWOFCli< 1-5) 

FNPFOR 

FOH 1 "-697. 726 
HARRIS! I >-l)X! 1-691 ) 

KNOT OR 

FOR I -727. 756 
HARRIS! I H-NFWOKIK 1-35 ) 

FNIMCIH 

RISK 

FOK 1*1.4 
HARRIS! 1 )-liX! 1 ) 

F.NllFOK 
FOK J-5. 8 

HARRIS! I )*NI W!iKIK 1-4 ) 

F.NllFOK 


o r 
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K)R 1*Y. 10 

Harris! i >*nx< i -a » 

ENtlSOR 

FOR I'll. 14 

HARRIS! 1 >=NEW0RIK 1-6 ) 

KNWON 

FOR 1 "1 M. 16 

HARRIS! 1 )-l IX< 1-8 ) 

FNUKOH 

FOR 1-17.04 

HARRIS! I l-NirWORO! 1-8 ) 

kNllKOR 

FOR 1 -OS. 07 

HARRIS! I 1-OX! 1-26 1 

F.NllKOR 

FOR 1-08. S27 

HARRIS! 1 I'-Nt.WORn! 1-11 ) 

fc N11FCR 

ENIHF 

IF! 4 1. t.Q. 1 1 TH£N 

1U-10S4 

EL. St 

IH'lMl 

tNlllF 

BUT-REROUT! ' 10. HARRIS. B. 1U. » 1ST AT. OWORJl ) 
CAli STATUS! '13 1 
t-NULCiOP 
FOR 1*1.0 

CALI XXYY! '13. '06 ) 

ENIffOK 

STOP 

ENU 


U i \ £ • 

<* POC«\v>, 


»w 

rr 


SU8R0UT 1 NK R I KNSLC RCOUNT , 1 N II . F 1 N. H. Nt- WORO ) 

C 

INTEGER* 6 HART 1 , R2. HO. B4, HR) . HH2. HHO. HH4. HH? 
ROUBLE HREC 1 " 1 ON S1C-N. F.XR, RANT , NKW0K1K 760 > 
i NT EGER H< 0000 ), SI ONI. FXP1 . INI I. KIN. NUMWRll. RCOUNT 
NtlMUHUM ! F 1 N- 1 N IT 1 J/4 
FOR 1-l.NUWUKH 

j.ini i>a*( i-i > 

SI ONI *H< J 1. AND. '700 
S 1 !»N 1 'S 1 C<N 1 SHIFT. -7 
IK! S 1 ON 1 KM. ) > 1HSN 
SION)— 1 
Kl SE 
SJONl'J 
K Nil IK 

SlUN-SIGNI 
HHJ'H! .M ) 

MH'^HRl AN) i ' 


700 




-38- 


Oi?- . . 

of pock g 


Al_l i V 


HH3*HH2. SHIFT. ! V 
ANJV '1 77 
BHb*F<H4 SHIFT. 2 6 
B2"KHb. OR. hhs 
H 3“H< J ♦'/ ). SHIFT, ft 
H4«R< J>3 > 

HANT 1 ■‘US'. CM?. W3. OR. B4 

MATO -MAN! 1 

IK MAN! 1. F.H. 0 ) THEN 

Nl-UOKIK KGOUNT HO. 

FXSfc 

fc.XPl«HI J ). AND. '177 
t.XHl«-fcXFl-VO 
fc.XP-16. ««E.XP1 

Nt.UOHIK RCXTUNT HS1GN*F.XP«HANT 
F.NII1F 

KCCHJN1 “KCCHJNT* 1 
F Niff- CUT 
RETURN 
34 b.Nl) 

C 

C: 

SUBROiniNfc ITRNSLC 1COUNT. IN IT. F I.*'. R. NEUVAR ) 

C 

INTEGFK«6 MANT1. B2. P3. P4. NfcWVARC 40 ). SION.*’, ftl. RBI PH2. RH3. Fft4, HI 
I NT EGER R< 3030 >. S I ON l . I NIT. PIN. NUF!WR1\ I COUNT 
NUPWHI1«< < FIN-IN1T HI >/4 
FOR 1-l.NUMV/Rl) 

J»1N1T>4«< 1-1 ) 

SIHXHBt J ). AND. '200 
SIGN1»S1GN1. SHIFT. -7 
1F< SIGN1. Ftt 1 ) THEN 

s;gni«^-i 

FL.SF. 

S1GN1-1 

F.NIJ1F 

V Sl(iN2«S]GNI 

RHHHCJ) p 

HHHH1 SHIFT. 2b 

RK2«H( J+ J ). ANIJ. ' 1 77 

HKHHC J+l ). AND. '200 

RR4-»RH2. SHIFT. 16 

HHb-HH3. SHIFT. 17 

B2* r HH4 CiN. HRi. 

H3"H< .H3 >. SHIFT, ft 
H4r-H< J+:i ) 

MANT HHl OR. » y OR. H3. OK K4 
NFWVAKC 1 (.01 IN I )=S1GN2*HANI 1 
) COUNT* ) CCll INT ♦ J 
FNIIFOK 
KF 1 1 H<N 
37 FNH 
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C 

SHHHMT i AKH S A 1 APH i‘KHATHH NY * .(HMTK 


(: 

ANH KHHUKMA IS ) 1 . 
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(: 

C 

CONTAIN IMP PAlA HkOM UNH HI IUHI TKACK 
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C; 
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C: 

HAT A HKfiHI > NASA 

PASSHS AKH CON'IAJNHP IN 


C 
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C: 
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C: 

(: 

C: 
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UOKLi 0 

HH SCRIPT ION 
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1 

HJUXJ 


C: 

> 

) PASS J 


C 

3 

MSHCX J ( ANT ) 


C 

r; 

4 

s 

MSHUXK AKC ) 
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c 
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C: 

7 

MJI.iXV 


(; 

8 

J PASSV 
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V 

MSHCX X ANT ) 


C 

10 
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(: 

1) 

KPX ANT I 


f; 

1> 

KPX AKC ) 


0 

JH 

t OH HAT A POINTS 


C: 
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1 AT 


C: 

784-SSH 

1 ON 


f: 

Si,4-B>3 

HI H.V 


C. 

874-10VH 

PHI H 


C: 
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CAI I STA1US< 'IV ) 

IK KIR I ). FU. 1 ) IHI'N 
SVh FOKMAK 1 X. H4. V. ) 
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x< y. < a ) 
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001 o soo 
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“i.l 

1*2 

U i»3 

"i,4 

"1.5 

"1.6 

1.7 

"i|8 



5 

7 

20 

1 

12 

19 

- 

- 

4 

3 

11 

20 

9 

45 

- 

- 

11 

21 

18 

20 

10 

41 

8 

- 

21 

35 

41 

45 

1 

39 

- 

- 

5 

10 

20 

48 

38 

19 

- 

6 

17 

19 

54 

44 

55 

39 

7 

20 

15 

13 

24 

47 

53 

35 

16 

19 

23 

30 

35 

61 

44 

36 

27 

30 

40 

46 

54 

54 

87 

86 

17 

20 

21 

32 

23 

39 

68 

109 

1 

- 

5 

10 

19 

18 

20 

67 

10 

10 

4 

3 

10 

12 

51 

88 

1 

10 

10 

8 

34 

59 

59 

54 

9 

- 

31 

53 

45 

38 

31 

66 

24 

48 

34 

34 

24 

25 

45 

56 

35 

37 

38 

29 

30 

42 

57 

43 

29 

17 

15 

19 

23 

40 

36 

41 

2 

13 

18 

10 

41 

45 

36 

78 

- 

k4 

6 

56 

35 

33 

29 

71 

10 

10 

43 

22 

42 

22 

12 

80 

— 

6 

27 

32 

14 

20 

26 

69 

4 

4 

- 

11 

14 

10 

71 

38 

— 

— 

— 

5 

14 

69 

40 

28 

— 

— 

— 

- 

10 

3 

45 

21 
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( 


0{ VGlXAL 

0F POOR 


PAG £ rs 

Quality 


1*9 

"i.lO 

1.11 

1.12 

1.13 

1.14 

1.15 

1,16 

27 

6 

10 

44 

26 

47 

65 

18 

31 

9 

6 

58 

24 

68 

46 

22 

55 

10 

- 

60 

44 

91 

10 

29 

38 

35 

1 

69 

62 

69 

20 

35 

41 

47 

39 

47 

100 

35 

31 

53 

47 

33 

78 

53 

80 

44 

64 

14 

83 

36 

77 

124 

40 

74 

67 

32 

64 

72 

104 

130 

77 

77 

61 

57 

68 

66 

183 

146 

129 

104 

107 

66 

92 

135 

311 

307 

275 

147 

73 

24 

113 

256 

306 

- 

- 

145 

155 

95 

171 

335 

- 

- 

- 

- 

152 

66 

115 

335 

- 

- 

- 

- 

109 

68 

77 

199 

159 

- 

- 

54 

76 

48 

108 

88 

136 

102 

53 

37 

37 

16 

99 

66 

88 

94 

20 

34 

37 

10 

106 

62 

53 

59 

71 

58 

38 

28 

51 

46 

51 

42 

49 

63 

51 

32 

44 

51 

35 

47 

19 

54 

44 

41 

42 

49 

14 

35 

40 

34 

38 

62 

46 

34 

1 

48 

40 

12 

39 

35 

43 

28 

- 

40 

40 

20 

23 

35 

28 

28 

5 

26 

32 

42 

2 

32 

27 

28 

10 

25 

23 

48 

2 

14 


( 




( 
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D 1.17 

“l.ie D 1.19 °1.20 D l,21 ®i,22 ®1.23 “l.24 

18 

1 9 4 9 10 


20 

14 

15 

21 

27 



_ 



29 

26 

31 

39 

4 





*■ 

35 

42 

21 

15 

23 

4 





46 

11 

29 

25 

15 

17 

3 




24 

28 

37 

31 

25 

19 

24 

16 



26 

46 

46 

38 

49 

45 

44 

8 



64 

68 

50 

24 

8 

14 

16 

17 



36 

21 

10 

15 

19 

17 

8 




30 

9 

18 

20 

21 

10 

9 

8 



72 

61 

54 

29 

17 

10 

10 

9 


i 

41 

51 

38 

48 

40 

49 

37 

22 


# 

* 

55 

21 

23 

29 

38 

21 

30 

32 


P 

32 

48 

22 

7 

6 

20 

21 

17 


■ 

28 

17 

35 

23 

12 

1 

3 

16 


: - 

33 

12 

18 

26 

23 

13 

6 



1 

21 

27 

6 

20 

21 

23 

10 

2 



— 

38 

25 

3 

22 

17 

21 

4 


I 

16 

19 

18 

17 

11 

33 

11 

10 



21 

14 

38 

11 

3 

2 

12 

4 


1 

1 

50 

17 

18 

20 

15 






31 

37 

19 

26 

13 

7 


_ 



28 

18 

27 

16 

11 

10 





34 

9 

5 


19 




* 

1 



a 
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TABLE II 


Summary of Characteristics of the "Unaffected" Passes 


PASS 

DATE 

TIME 

KP 

AVG ELEV 

49 

11/ 5/79 

5:45 

0 

573 

50 

11/ 5/79 

7:19 

3 

577 

51 

11/ 5/79 

8:53 

3 

578 

52 

11/ 5/79 

10:27 

7 

578 

75 

11/ 6/79 

22:24 

3 

577 

76 

11/ 6/79 

23:58 

3 

576 

200 

11/15/79 

1:42 

0 

573 

260 

11/18/79 

23:21 

3 

562 

262 

11/19/79 

2:29 

3 

562 

263 

11/19/79 

4:02 

7 

560 

308 

11/22/79 

2:15 

0 

551 

309 

11/22/79 

3:49 

0 

550 

323 

11/23/79 

1:39 

3 

547 

325 

11/23/79 

4:46 

0 

546 

400 

11/28/79 

1:44 

3 

522 

401 

11/28/79 

3:18 

0 

522 

402 

11/28/79 

4:51 

0 

521 

411 

11/28/79 

18:53 

7 

518 

413 

11/28/79 

22:00 

3 

517 

414 

11/28/79 

23:34 

3 

517 

415 

11/29/79 

1:07 

3 

517 

416 

11/29/79 

2:41 

3 

516 

417 

11/29/79 

4:41 

0 

516 

418 

11/29/79 

5:48 

0 

515 

419 

11/29/79 

7:21 

7 

515 

420 

11/29/79 

8:55 

7 

515 

757 

12/21/79 

5:40 

0 

396 

821 

12/25/79 

9:10 

0 

387 

1019 

1/ 7/80 

4:48 

0 

385 

1168 

1/16/80 

20:02 

3 

415 

1169 

1/16/80 

21:35 

3 

415 

1251 

1/22/80 

4:45 

0 

438 

1282 

1/24/80 

4:48 

0 

448 

1298 

1/25/80 

5:36 

0 

452 

1299 

1/25/80 

7:09 

7 

453 

1300 

1/25/80 

8:42 

7 

453 

1461 

2/ 4/80 

18:08 

3 

494 

1463 

2/ 4/80 

21:14 

3 

492 

1464 

2/ 4/80 

22:47 

- 3 

495 

1465 

2/ 4/80 

0:20 

0 

495 

1466 

2/ 5/80 

1:53 

0 

495 

1467 

2/ 5/80 

3:26 

0 

495 

1468 

2/ 5/80 

4:59 

0 

496 
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1469 

2/ 5/80 

6 

32 

7 

496 

1556 

2/10/80 

21 

12 

7 

507 

1557 

2/10/80 

22 

45 

7 

507 

1559 

2/11/80 

1 

50 

3 

508 

1575 

2/12/80 

2 

36 

7 

509 

1576 

2/12/80 

4 

08 

0 

509 

1592 

2/13/80 

4 

54 

0 

510 

1593 

2/13/80 

6 

26 

0 

510 

1594 

2/13/80 

7 

59 

0 

510 

1695 

2/19/80 

20 

10 

3 

506 

1775 

2/24/80 

23 

48 

3 

496 

1868 

3/ 1/80 

23 

24 

3 

475 

1880 

3/ 2/80 

17 

56 

7 

472 

1889 

3/ 3/80 

7 

49 

3 

470 

1903 

3/ 4/80 

5 

25 

3 

466 

1928 

3/ 5/80 

19 

59 

7 

459 

1929 

3/ 5/80 

21 

32 

7 

458 

1930 

3/ 5/80 

23 

04 

7 

456 

2008 

3/10/80 

23 

20 

7 

434 

2014 

3/11/80 

8 

35 

3 

432 

2025 

3/12/80 

1 

32 

0 

429 

2026 

3/12/80 

3 

05 

0 

428 

2027 

3/12/80 

4 

37 

0 

428 

2028 

3/12/80 

7 

42 

0 

427 

2030 

3/12/80 

9 

14 

0 

427 

2033 

3/12/80 

13 

52 

0 

426 

2034 

3/12/80 

15 

24 

7 

426 

2035 

3/12/80 

16 

57 

7 

425 

2073 

3/15/80 

3 

29 

0 

414 

2074 

3/15/80 

5 

02 

0 

413 

2075 

3/15/80 

6 

34 

3 

413 

2076 

3/15/80 

8 

06 

3 

413 

2079 

3/15/80 

12 

44 

3 

412 

2080 

3/15/80 

14 

16 

3 

412 

2081 

3/15/80 

15 

48 

3 

4 il 

2082 

3/15/80 

17 

21 

3 

411 

2083 

3/15/80 

18 

53 

0 

411 

2084 

3/15/80 

20 

26 

0 

411 

2134 

3/19/80 

1 

25 

3 

396 

2135 

3/19/80 

2 

57 

3 

396 

2213 

3/24/80 

2 

59 

3 

376 

2214 

3/24/80 

4 

31 

7 

375 

2230 

3/25/80 

5 

07 

0 

372 

2231 

3/25/80 

6 

40 

3 

372 

2232 

3/25/80 - 

8 

12 

3 

371 
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ORIGINAS- . . ... > 

OF POOR QUALiTY 



-5 

-4 


-4 

-2 

-2 

1 

0 

-1 

-2 

-2 

-2 

-2 

1 

-1 

1 


1 


TABLE III 


24 z 24 Array AB. . of Scalar Anonollaa 

*ij 








-2 

-3 

-3 

-3 

-3 

-2 

-2 

1 

0 

-1 

-1 

1 

0 

0 

2 

2 

1 

1 


-6 

-2 

-2 

-2 

1 

0 

-1 

0 

-1 

1 

0 

1 

2 

1 

1 

0 

-1 

1 

2 

1 

1 


-5 

-1 

-1 

-2 

-2 

1 

3 

4 
2 
2 
0 
1 
2 
3 
2 
0 

-1 

-1 

1 

1 

0 

-1 


2 

0 

-1 

-2 

1 

1 

2 

3 

4 
1 
1 
1 
2 
2 
2 
0 

-2 

-2 

0 

1 

-1 

-1 

0 

1 


14 

2 

-5 

-1 

1 

1 

1 

2 

3 
2 

-1 

4 
4 
2 
0 
1 

-1 

-2 

0 

0 

0 

1 

2 

2 


7 

4 

-3 

-3 

-1 

-1 

2 

4 

3 

2 

2 

3 

3 

3 

3 

1 

-1 

-2 

-2 

1 

0 

1 

0 

-2 



6 

1 

-1 

-3 

-1 

0 

0 

5 

4 

1 

1 

-1 

-1 

3 

4 
2 

-1 

-2 

0 

1 

1 

1 

1 

-1 



ft 


<5 


ORIGINAL PAtU 4 
OF PC OR QUALITY 



(9 
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ORIGINAL •: 
OF POOR QUALITY 


AB 


AB 


AB 


AB. 


AB. 


1,17 1,18 1,19 1,20 1,21 


AB 


1.22 


0 

0 

4 

2 

-2 

2 

3 

4 
4 
3 

-1 

-2 

1 

3 

4 
8 
4 

-2 

-7 

-2 

4 

4 

4 


4 
2 
3 
1 

-1 

1 

6 

6 

5 
0 

3 

4 
2 

5 
7 

7 

8 
5 

-6 

-6 

0 

3 

4 

5 


-2 

1 

3 
2 

-1 

1 

4 

4 

5 
-3 

5 

5 

-3 

4 

5 
8 
0 

-5 

-6 

-5 

-3 

6 
5 
2 



-4 

-3 

-7 

-5 

-3 

-5 

-2 

-3 

-5 

-3 

1 

-1 

3 

2 

1 

-3 

-4 

-6 

-3 

-7 

-3 

2 

0 

1 


-3 


-2 

-2 

-4 

0 

-1 

-3 

2 

2 

2 

-2 

1 

-3 

-3 

-3 

-4 

-3 


2 

5 


AB 


1,23 ^* 1,24 


-4 

-1 

1 

-1 

-2 

2 

0 

-1 

-4 

-2 

-1 

-3 

-2 

-3 

-1 

1 


-5 

-1 

0 

1 

3 

-2 

-2 

-5 

-5 

0 


-2 

4 

6 



ORIGAMI t'AuL i ,j 

OF POOR QUALITY 


TABLE IV 

24 z 24 Array AX. of Radial Anoaollaa 
lfj 



i 


t 


ORIGINAL 

OF POOF 
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ORIGINAL P/SQ£ ts 
OF POOR QUALITY 


AX 

i .17 

AX 

1,18 

AX 

1,19 

**1 

6 

27 

2 

2 

7 

2 

8 

1 

0 

5 

0 

-2 

-2 

-2 

-1 

-1 

-6 

-6 

-8 

1 

-13 

-14 

2 

0 

8 

24 

7 

6 

10 

6 

21 

16 

0 

-13 

-2 

13 

-14 

-13 

-9 

-34 

-13 

-14 

-7 

4 

5 

-8 

-14 

-10 

-15 

-12 

-5 

-20 

-3 

7 

-29 

-40 

-10 

4 

3 

-27 

6 

-13 

6 

-12 

-1 

-10 

-8 

4 

- 

-9 

4 

-4 

7 

-8 

14 

-2 

-3 

-26 

11 

40 

5 

-39 

11 

-5 

-11 

-10 

4 

3 

-4 

3 

17 

-6 

-1 

3 

1 

- 


AX 

i .21 

AX 

1,22 

AX 

1.23 

“ i . 

6 

-1 

m 

. 

1 

— 

m 

m i 

-10 

• 

mm 

m 

2 

17 

m 

m 

-2 

-12 

-28 

_ 

-8 

-8 

-2 

-3 

2 

-5 

-14 

-24 

3 

1 

-6 

-7 

26 

3 

-5 

— 

-3 

29 

22 

23 

0 

-5 

34 

9 

-7 

1 

11 

7 

-8 

23 

28 

17 

9 

17 

23 

-2 

31 

-19 

51 

8 

1 

18 

12 

— 

-17 

27 

10 

4 

5 

44 

3 

-7 

-21 

-4 

-5 

-5 

17 

4 

1 

-11 

-2 

2 

- 

_ 

-6 

-8 

- 

_ 

4 

-8 

- 

- 

-4 

- 
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ORIGINAj. 'V.CE fS 
OF POOR QUAUTY 


TABLE V 

24 z 24 Array AY. . of TangontUl Anoxalio* 
*»j 



i 


a 


ORIGINAL P/.^t: 13 

OF POOR QUALITY 



d 


OF 



AY 1,18 

AI 1,19 

" i . 

-3 

-5 

2 

-2 

0 

2 

1 

-1 

-2 

-3 

-1 

1 

-7 

-6 

2 

-3 

-8 

-6 

-11 

-8 

-18 

-18 

4 

-1 

7 

12 

7 

-1 

5 

-4 

7 

-4 

- 

5 

12 

- 

-12 

-18 

15 

7 

3 

1 

3 

7 

-3 

-2 

9 

5 

19 

11 

17 

9 

13 

9 

1 

6 

9 

7 

5 

3 

-6 

13 

5 

22 

0 

3 

20 

-11 

- 

16 

26 

15 

26 

37 

-6 

-1 

10 

-14 

-13 

-21 

-13 

-12 

-15 

-12 

-7 

-17 

1 

-6 

-7 

-6 


-8 

-9 

-18 


- 
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0 
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i t * % 
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TABLE VI 


24 z 24 Array of A2 4 of Vertical Anomalies 
*»j 


^i.l “i.2 ^i.3 ^i.4 ^i.5 ^i.6 ^i.7 ^i.8 


8 

5 

4 

4 

4 
1 
1 
1 

3 
2 

5 

4 
1 
4 

0 

0 


5 

2 

2 

2 

4 

5 
5 
0 

1 

0 

2 

0 

1 

1 

4 

-2 

-3 

0 


9 

5 

1 

4 

0 

2 

4 

2 

1 

0 

0 

-l 

0 

-2 

-1 

-1 

1 

-1 

-1 

1 

0 

0 


11 

3 

5 

5 

1 

0 

-1 

-1 

-2 

-3 

0 

1 

-5 

-3 

0 

1 

0 

-1 

1 

3 

3 


2 

-1 

0 

2 

0 

-1 

-3 

-2 

-5 

0 

1 

-4 

-4 

-4 

-L 

2 

1 

-1 

0 

3 

5 

3 

-2 


12 

-3 

3 

2 

-1 

-3 

-3 

-3 

-3 

-3 

1 

-3 

-5 

-5 

-3 

-3 

-2 

-1 

-1 

-1 

0 

2 

-1 

0 


-3 

-5 

4 

4 

2 

-2 

-3 

-4 

-4 

-4 

-2 

-5 

-5 

**6 

-5 

-4 

-3 

0 

-1 

-1 

-1 

-1 

0 

1 


-3 

-1 

-1 

3 

2 

0 

-2 

-5 

-6 

-3 

0 

-2 

-2 

-5 

-b 

'5 

-1 

2 

0 

-2 

-2 

-2 

-1 

2 
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0RIG1NAL P. - iS 
OF POOR QUALITY 




hi 


1.10 “l.U “1.12 “1.13 “1.14 ”1.15 “1.16 


hi 


hi 


hi. 


hi . 


A2 J 


0 

2 

4 

4 

1 

1 

2 

2 

2 

1 

0 

3 

1 

3 

3 

4 
1 
2 
1 
1 
2 
2 
0 
1 


-1 

2 

5 

5 
1 

-2 

-5 

-4 

-2 

0 

6 
7 
3 

-3 

-1 

-3 

2 

1 

1 

-1 

-4 

-2 

1 

5 


2 

5 

-4 

-3 

-2 

-1 

1 

0 

-1 

7 


-1 

0 

3 

5 

5 

6 

-1 

-4 

2 

1 


-2 

3 

4 

-2 

-4 

-3 

1 

1 

-2 

-2 


1 

4 

4 

-2 

-3 

-5 

-6 

-4 

-1 

2 


0 

5 

1 

-7 

-7 

0 

4 

2 

0 

-3 


0 

2 

1 

2 

1 

0 

-5 

-2 

1 

4 


4 

4 

-1 

-6 

-4 

0 

2 

1 

2 

0 

10 


-8 

-1 

2 

0 

2 

5 

2 



-4 

0 

4 


3 
1 

-2 

-3 

-2 

-2 

-3 

-1 

2 

-1 

4 
6 
1 

-2 

-1 

c 

1 

1 

3 

8 

3 

-2 

1 

8 


2 

3 

-1 

-2 

0 

-2 

-7 

-7 
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-1 

5 

7 

3 

-4 

-5 

-1 

-5 

0 

3 

3 

1 

-4 

-2 

5 
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